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Abstract 

A scheme for making ab-initio calculations of the dynamic paramagnetic spin 
susceptibilities of solids at finite temperatures is described. It is based on Time- 
Dependent Density Functional Theory and employs an electronic multiple scattering 
formalism. Incommensurate and commensurate anti-ferromagnetic spin fluctuations 
in paramagnetic Cr and compositionally disordered CV95V5 and Crg^Re^ alloys are 
studied together with the connection with the nesting of their Fermi surfaces. We 
find that the spin fluctuations can be described rather simply in terms of an over- 
damped oscillator model. Good agreement with inelastic neutron scattering data is 
obtained. 
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Chromium is the archetypal itinerant anti-ferromagnet (AF) whose famous incommensurate 
spin density wave (SDW) ground state is determined by the nesting wave- vectors q ne st identified in 
the Fermi surface M. Chromium alloys also have varied AF properties |1] and their paramagnetic 
states have recently attracted attention owing, in part, to analogies drawn with the high temperature 
superconducting cuprates especially (La c SVi_ c ) 2 Cw04. The incommensurate SDW fluctuations in 
these materials are rather similar to those seen in the paramagnetic phase of Cr close to T^. 
Moreover 'parent' materials La^CuO^ in the one instance and Cr^Mn^ or Cr^Re^ in the other 
are simple commensurate AF materials which on lowering the electron concentration by suitable 
doping develop incommensurate spin fluctuations which may be promoted by imperfectly nested 
Fermi surfaces. 

Here we examine the nature of damped diffusive spin fluctuations in chromium above the Neel 
temperature T/v which are precursory to the SDW ground state. We also study dilute chromium 
alloys, Cr 95 i?e5 and Cr 95 V^ and obtain good agreement with experimental data. For example, recent 
inelastic neutron scattering experiments have measured incommensurate AF 'paramagnons', 
persisting up to high frequencies in the latter system. We explore the temperature dependence, 
variation with dopant concentration and the evolution of the spin fluctuations in these systems from 
incommensurability to commensurability with increasing frequency and provide the first ab-initio 
description of these effects. To this end we describe a new scheme for calculating the wave- vector and 
frequency dependent dynamic spin susceptibility of metals which is based on the Time Dependent 
Density Functional Theory (TDDFT) of Gross et al. |4| and as such is an all electron theory. For the 
first time the temperature dependent dynamic spin susceptibility of metals and alloys is calculated 
from this basis. There have been several simple parameterised models to describe the magnetic 
properties of Cr and its alloys f7j. All of these have concentrated on the approximately nested 
electron 'jack' and slightly larger octahedral hole pieces of the Fermi surface U and, at best, have 
only included the effects of all the remaining electrons via an electron reservoir. We find similarities 
between our results and results from such models but show that a complete picture is obtained only 
when an electronic band-filling effect which favors a simple AF ordering at low temperature is also 
considered. We also find that the spin fluctuations are given an accurate description as overdamped 
diffusive simple harmonic oscillator modes which are at the heart of theories of the effects of spin 



fluctuations upon the properties of itinerant electron systems [Q. 

Over the past few years great progress has been made in establishing TDDFT [|J. Analogs 
of the Hohenberg-Kohn || theorem of the static density functional formalism have been proved 
and rigorous properties found. Here we consider a paramagnetic metal subjected to a small, time- 
dependent external magnetic field, b(r, t) which induces a magnetisation m(r, t) and use TDDFT 
to derive an expression for the dynamic paramagnetic spin susceptibility x(q> w ) v i & a variational 
linear response approach ||. Accurate calculations of dynamic susceptibilities from this basis are 
scarce (e.g. 0) because they are difficult and computationally demanding. Here we mitigate these 
problems by accessing x(q, w) via the corresponding temperature susceptibility x(q, w n ) where w n 
denotes a bosonic Matsubara frequency flQfl . We outline this approach below. 

The equilibrium state of a paramagnetic metal, described by standard DFT, has density Po( r ) 
and its magnetic response function x(rt; r't') = (5m[b}(r, t)/Sb(r', t'))|ft =0iPo is given by the following 
Dyson-type equation. 

X(rt;r't') = X s(rt;r't') + J <fn J dh J dr 2 J dt 2Xs (rt; r^K^ih; T 2 t*)x(T % t 2 , r't') (1) 

Xs is the magnetic response function of the Kohn-Sham non-interacting system with the 
same unperturbed density po as the full interacting electron system, and K xc (rt; r't') = 
(5b xc (r,t)/5m(r',t'))\f, =0t p is the functional derivative of the effective exchange-correlation mag- 
netic field with respect to the induced magnetisation. As emphasised in ref. [§] eq.l represents an 
exact representation of the linear magnetic response. The corresponding development for systems 
at finite temperature in thermal equilibrium has also been described . In practice approximations 
to K xc must be made and this work employs the adiabatic local approximation (ALDA) jp so that 
Ki LDA {vt;r't') = (rf 2 6^(p(r,t),m(r,t))/rfm 2 (r,t))| po , m=0 5(r-r')5(t-t') = I xc (r)5(r-r')5(t-t'). 
On taking the Fourier transform with respect to time we obtain the dynamic spin susceptibility 
X(r,r';w). 

For computational expediency we consider the corresponding temperature susceptibility [|TD[ 
x(r, r'; w n ) which occurs in the Fourier representation of the temperature function x( rr ; r ' r ') that 
depends on imaginary time variables t,t' and w n are the bosonic Matsubara frequencies w n = 
InnksT . Now x(r,r';w n ) = x(r, r'; iw n ) and an analytical continuation to the upper side of the 



real w axis produces the dynamic susceptibility x(r, r'; w). Using crystal symmetry and carrying out 
a lattice Fourier transform we obtain the following Dyson equation for the temperature susceptibility 

X(x, x, q, w n ) = x s (x, x', q, w n ) + J dxix s (x, xi, q, w n )/ :EC (xi)x(xi, x', q, w n ) (2) 

with x,x' and xi measured relative to crystal lattice unit cells of volume Vws- 

In terms of the DFT Kohn-Sham Green function of the static unperturbed system 

% s (x, x', q, w n ) = —afiTr. ^ G ( x ' x '> R > f 1 + iv m)G(x!, x, -R, fi + %{v m + w n ))e iq ' R (3) 

P R m 

where R is a lattice vector between the cells from which x and x' are measured, /i the chemical 
potential and fermionic Matsubara frequency (2n + l)nkBT. The Green function can 

be obtained within the framework of multiple scattering (KKR) theory [11]. This makes this 
formalism applicable to disordered alloys as well as ordered compounds and elemental metals, the 
disorder being treated by the Coherent Potential Approximation (CPA) ||12|| . Then the partially 
averaged Green function, (G(r, r', z)) ra ,r'i3, where r,r' lie within unit cells occupied by a and (3 
atoms respectively, can be evaluated in terms of deviations from the Green function of an electron 
propagating through a lattice of identical potentials determined by the CPA ansatz [|13|]. 

To solve equation (|2]), we use the direct method of matrix inversion and local field effects are 
fully incorporated. x(q, q; w n ) = (I /Vws) I d x I rfx'e* q '^ x_x '^x( x ? x '? w n) can then be constructed. 
The most computationally demanding parts of the calculation are the convolution integrals over 
the Brillouin Zone which result from the expression for Xs, eq. @. Since all electronic structure 
quantities are evaluated at complex energies, these convolution integrals have no sharp structure 
and can be evaluated straightforwardly by an application of adaptive quadrature [14[ . 



As discussed in ref. |I(J, for example, we can define the retarded response function x(q, q, z) 
of a complex variable z. Since it can be shown JTO] formally that lim^oo x( z ) ~ ^/z 1 and we 
can obtain x(zw n ) from the above analysis it is possible to continue analytically to values of z 
just above the real axis, i.e. z = w + irj. In order to achieve this we fit our data to a rational 
function x(q, q,w n ) = x(q)(l + Lfcl"i 2 Uk(q}w^)/(1 + J2k=i -^fc(q)w^) with the choice of coefficients 
Uk,Dk ensuring that the sum rule involving the static susceptibility x(q) is satisfied, i.e. x(q) = 
(2/tt) Jq°° dwlmx(q., q, w)/w. We find that very good fits are obtained with small M. 



For chromium and its alloys, we find that M = 2 is sufficient to provide excellent fits to the 
calculations of x over a wide range of w n 's, i.e. X _1 (q, q, w n ) = x~ 1 (<i)(l + (w n /r(q)) + (w n /Q(q)) 2 ) 
so that x _1 (q, q,w) = x _1 (q)(l — i(u>/T(q)) — (w/Q(q)) 2 ) (standard error < 3% of mean). For 
the systems studied here we find f2(q)/T(q) < 0.15 <C 1 and so the spin dynamics can be described 
in terms of a heavily overdamped oscillator model. Evidently £s.p(q) = h/T(q) represents a relax- 
ation time for a damped diffusive spin fluctuation with wavevector q. Moreover, the imaginary 
part of the dynamical susceptibility which, when multiplied by (1 — exp(— ftw))' 1 , is proportional 
to the scattering cross-sections measured in inelastic neutron scattering experiments, is written 
Jmx(q, q, w) = x(q)wr _1 (q)/((l — (w/Q(q)) 2 ) 2 + (u>/T(q)) 2 ). We note that theories for the spin 
fluctuation effects upon itinerant electron properties, including quantum critical phenomena, also 
invoke such a model 0. The small Q, = (q — q ), dependence of X _1 (Qq + Q) an d r(q + Q) is of 
particular importance. (q is where x -1 (q) is smallest.) 

Finite-temperature calculations were carried out for the static susceptibilities of the three sys- 
tems, using the experimental b.c.c. lattice spacing of Cr, 2.88 A, and von Barth-Hedin local 
exchange and correlation ||15[ . We find that (i) Cr orders into an incommensurate AF state below 
280K specified by q = {0,0,0.93}, where experiment yields T N = 311K and q = {0,0,0.95} [|TJ; 

(ii) CV95V5 does not develop magnetic order at any temperature, as found in experiment [|l]|§; an d 

(iii) Cr 95 Re 5 orders into a weakly incommensurate AF state below T=410K (q = {0,0,0.96}), 
whereas experimentally it forms a commensurate AF state below TV of 570K 

Our calculations for 7mx(q, q, w) are shown in figures 1(a) and (b) for CV95V5 and Cr$§Re§ 
respectively. Our calculation of 7mx(q, q, w) for paramagnetic Cr at T = 300-fT is broadly similar 
to that for paramagnetic Cr at T = OK by Savrasov || so a figure is not presented. It shows 
incommensurate spin fluctuations for small frequencies which are signified by peaks in Jmx(q, q, w) 
at q which is equal to the Fermi surface nesting vector q nest . For increasing w the peaks move to 
q = {0, 0, 1} i.e. the spin fluctuations become commensurate. The spin fluctuations at 300K shown 
in fig. 1(a) for CV95V5, on the other hand, remain incommensurate up to much higher frequencies 
maintaining intensity comparable to that at the peak at low w. This qualitative difference between 
the two systems has not been described before by a first-principles theory although found exper- 
imentally For lower temperatures we find that Imx(q, q, w) of Cr 95 V 5 becomes a sharper 



function of w and we can also infer that (1 — exp(— (3w))~ l Imx{Q., q, w) should vanish for small w 
when T — > OK. These aspects have also been noted from experimental measurements [|lj. 

It is striking that the alloy's Fermi surface is well defined [ETIJ despite impurity scattering although 
it is more poorly nested (the difference between the sizes of the electron and hole octahedra is larger) 
than that of Cr owing to its fewer electrons. Once again the peaks in 7mx(q, q, w) occur at the 
nesting vectors q nes t = {0,0,0.9}. The spin fluctuations in the paramagnetic phase of Cr 95 Re 5 
are shown in fig. 1(b). Here adding electrons by doping with Re has improved the Fermi surface 
nesting so that q nes t = {0, 0, 0.96} and imx(q, q, w) has weight spread from q nest to {0, 0, 1}. The 
dominant spin fluctuations now rapidly become commensurate with increasing w. We obtain a 
rather similar picture from the calculations for Cr by artificially raising the chemical potential by a 
small amount. Interestingly when we account for thermally induced electron-phonon scattering 
by adding a small shift (« 20 meV) to the Matsubara frequencies v m in eq. (|3]), we find a tendency 
for the dominant spin fluctuations to become commensurate at lower w in both Cr and Cr^Re^. 

Some of these features also emerge qualitatively from the simple parameterised models based on 
that part of the band-structure near \i which leads to the nested electron and hole octahedra at the 
Fermi surface 0. Our 'first-principles' calculations, being based on an all-electron theory, however, 
need some additional interpretation. As analysed by recent total energy calculations [|l(J, b.c.c. Cr 
with the experimentally measured lattice spacing tends to form a commensurate AF phase at low 
temperatures which is modulated by a spin density wave of appropriate wavelength. The overall 
AF instability of the paramagnetic phase is promoted by the approximate half- filling of the narrow 
3d-bands [[H]] which is further modified by a weak perturbation coming from the Fermi surface 
nesting. As dopants such as V and Re are added not only is the Fermi surface nesting altered but 
also the d-bands become either further from or closer to being half-filled. 

The calculations can be summarised in terms of the damped oscillator model. x _1 (qo + Q) — 
X _1 (qo + cQ 2 ) for small Q for Cr and CY95V5 whereas for Crg 5 Re$, X _1 (q) is nearly constant for a 
range of q between q nest and {0, 0, 1}. We find the product 7(q) of x(q) with damping factor T(q) 
to be only very weakly temperature dependent in these three systems and ~ 7(qo) ; a constant, 
for small Q, yielding a dynamical critical exponent || of 2 typically assumed for antiferromagnetic 
itinerant electron systems. The nature of the spin fluctuations can be succinctly described via the 



variance < m 2 (q) >. From the fluctuation dissipation theorem, < m 2 (q) >= (1/vr) j^dwil — 
exp(— f3w))^ 1 Imx((l, q, w). Fig. 2 shows < m 2 (q) > at several temperatures for Cr where we have 
used a frequency cutoff of 500 meV and so have not included the faster of the quantum fluctuations. 
Near T/v the magnetic fluctuations have their greatest weight around the q nes t- At higher T the peak 
diminishes and weight grows at q's nearer {0,0, 1} reflecting the shift in the peak in 7mx(q, q, w) 
from (\ nes t to commensurate q's with increase in frequency w. Similar plots to fig. 2 for CV95V5 and 
CrgsRes show respectively a smaller and greater tendency for the weight in < m 2 (q) > to transfer 
in this way. If the frequency cutoff is reduced, < m 2 (q) > near {0,0, 1} is sharply diminished so 
that the Brillouin zone integral of < m 2 (q) >, < m 2 >, decreases with increasing temperature as 
inferred from neutron scattering data [l|]. 

We have not included the effects of spin fluctuation interactions i.e. mode-mode coupling |18| 



into our calculations and have determined T/v and the static susceptibility by what is essentially an 
ab-initio Stoner theory. In weak itinerant ferromagnets, for example, mode-mode coupling causes 
a dramatic suppression of the Curie temperatures from those estimated from a Stoner theory. 
For the Cr systems studied here, the fair agreement with experiment which we obtain for T/v 
and the relatively large value of the damping factor T(q) with respect to that in weakly itinerant 
ferromagnets, is suggestive that such spin fluctuation effects are small. Spin fluctuation calculations 
have, however, been carried out by Hasegawa and others for simple parameterised models of Cr 
neglecting Stoner particle- hole excitations and Fermi surface nesting Jl9| . Using a functional integral 
technique he made a high temperature (static) approximation so that all the thermally induced 
fluctuations were treated classically and found T/v for commensurate AF order to be 370K and 



V < m 2 > to increase linearly with temperature above T/v- A quantitative calculation, however, 
in which the Stoner particle-hole excitations and spin fluctuations are treated within the same 
framework is needed to determine unequivocally whether or not a Stoner-like picture is adequate 
for these systems. 

In summary, we have presented a first-principles framework for the calculation of dynamic 
paramagnetic spin susceptibilities of metals and their alloys at finite temperatures. At this point 
we add that the approach can also be applied to the study of magnetic excitations in magnetically 
ordered materials. The first applications on the AF spin fluctuations in Cr and Cr 95 Re 5 above T N 



and in paramagnetic Cr 95 V 5 have found good agreement with available experimental data. 
We are grateful to F.J.Pinski, S.Hayden and R.Doubble for useful discussions. 
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FIGURES 

FIG. 1. /mx(q,q,«)) of (a) Cr 95 V 5 at T =300K and (b) Cr 95 Re 5 at T = 1MT N (450K) in units of 
fj^eV" 1 for wave- vectors q along the {0, 0, 1} direction where q is in units of it /a (a is the lattice spacing). 
The frequency axis, w, is marked in meV. 

FIG. 2. The variance of the spin fluctuations in Cr < m 2 (q) > in jj? B for wave- vectors q along the 
{0,0, 1} direction at 350K (full line),400K (long dashes), 450K (medium dashes) and 700K (short dashes). 
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